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The six Higgs doublet model is a minimal extension of the Standard Model (SM) that addresses 
dark matter and gauge coupling unification. Another Higgs doublet in the 5 representation of a 
discrete symmetry group, such as Se, is added to the SM. The lightest components of the 5-Higgs 
are neutral, stable and serve as dark matter so long as the discrete symmetry is not broken. Direct 
and indirect detection signals, as well as collider signatures are discussed. The five-fold multiplicity 
of the dark matter decreases its mass and typically helps make the dark matter more visible in 
upcoming experiments. 



I. INTRODUCTION 

The Standard Model (SM) of particle physics has 
enjoyed great success in explaining physics below the 
electroweak scale, but it is unlikely to remain the sole 
description of nature up to the Planck scale. The SM 
does not contain a viable cold dark matter candidate, 
such as a weakly interacting massive particle (WIMP), 
or satisfactorily address the issue of gauge coupling 
unification. Both of these issues hint at new physics 
beyond the electroweak scale. In addition, the SM must 
be fine-tuned to regulate the large radiative corrections 
to the Higgs mass and the cosmological constant. It is 
possible to obtain a natural Higgs from supersymmetry 
PP, large extra dimensions |2], technicolor [3J, Randall- 
Sundrum models |4], or the little Higgs mechanism [5]. 
However, these models do not address the magnitude 
of the cosmological constant, which appears to be more 
fine-tuned than the Higgs mass, and leads us to question 
the central role of naturalness in motivating theories 
beyond the Standard Model. 

An alternative approach is to explore non-natural 
extensions of the SM in which fine-tuning is explained 
by environmental selection criteria [BJ. Weinberg noted 
that if the cosmological constant was much larger than 
its observed value, galaxy formation could not occur 
[7]. Additionally, the Higgs vacuum expectation value 
cannot vary by more than a factor of a few before atoms 
become unstable [HI IS]- In the context of the string 
theory landscape populated by eternal inflation, there 
exists a natural setting for environmental selection to 
play out [HI Q2]. This motivates considering models 
with parameters that may not be natural, but which are 
forced to be small by environmental selection pressure. 

Split supersymmetry is an example of a model that 
relaxes naturalness as a guiding principle and focuses 
on unification and dark matter [9j [lOj HTj - In this 
theory, the scalar superpartners are ultra-heavy and the 



electroweak scale consists of one fine-tuned Higgs and 
the fermionic superpartners, which are kept light by R- 
symmetry. The additional fermions alter the running of 
the gauge couplings so that unification occurs at high 
scales. 

The "minimal model" was presented in O ITS] , 
where a Dirac electroweak doublet serves as a dark 
matter candidate and leads to gauge coupling unifica- 
tion. A fermion singlet that mixes with the dark matter 
must also be introduced to avoid conflicts with direct 
detection results. The additional fermion leads to richer 
phenomenology, but at the cost of introducing a new 
mass scale. 

In both these models, fermions serve as dark matter 
and technical naturalness protects their masses from 
large radiative corrections. Split supersymmetry as- 
sumes that the high energy dynamics are supersym- 
metric, but that high-scale susy breaking is preferred. 
If the susy breaking sector communicates R-symmetry 
breaking inefficiently, the gauginos and Higgsinos end up 
much lighter than the typical supersymmetric particles. 
This should be contrasted with models like those in [14 , 
where an ad hoc dynamical mechanism is invoked to 
make the fermionic dark matter much lighter than the 
GUT scale. Without understanding how the "minimal 
model" fits into a high energy theory, it may be that it 
requires fine-tunings of fermion masses to get a viable 
dark matter particle. Similarly, without understanding 
how R-symmetry breaking is explicitly communicated to 
the gauginos and Higgsinos, it may be that split susy 
requires a tuning of a fermion mass to get weak-scale 
dark matter. 

The use of technical naturalness to justify new light 
fermions may be particularly misleading for relevant 
couplings that determine the large-scale structure of the 
Universe. In [TBJ[T7], the formation of galactic structures 
with properties similar to the Milky Way places bounds 
on the ratio of the dark matter density to baryon density. 
Typically, these bounds are not as strong as those on the 



cosmological constant, but they fix the dark matter mass 
to within an order of magnitude [16] (or three orders of 
magnitude in 17J). This opens up the possibility that 
the mass of the dark matter is unnatural and is set by 
environmental conditions so that the baryonic fraction 
of matter is not over- or under-diluted. 

A candidate for unnatural dark matter is the scalar 
WIMP. While there has been some recent work on 
minimal models with scalar dark matter [T3 H3 EE I5U1 
\21\ 122] . none provide a framework for gauge unification. 
In this paper, we will study a minimal scalar unifon 
sector that also contains a viable dark matter candidate. 
An additional Higgs doublet is added to the SM that is 
in a 5 or 6-plet of a new global discrete symmetry* This 
global symmetry remains unbroken, yielding a spectrum 
of two five-plets of real neutral scalars and one five- 
plet of charged scalars. Relic abundance calculations 
give the WIMP mass to either be ~ 80 GeV or in the 
range ~ 200-700 GeV. The model will be tested by next- 
generation direct and indirect detection experiments, 
and may possibly have signatures at the LHC. 

The model will be presented in greater detail in Sec. 
II. In Sec. Ill, the renormalization group equations are 
solved to illustrate gauge unification and the allowed 
weak-scale values of the theory's quartic couplings. The 
relic abundance calculation is presented in Sec. IV 
and the predicted experimental signals, in Sec. V. The 
results are summarized in Sec. VI. 



II. THE SIX HIGGS DOUBLET MODEL 

The model proposed in this paper consists of the 
Standard Model Higgs, h, plus an additional electroweak 
doublet, i?5, in the 5 representation of a global discrete 
symmetry group. The discrete symmetry is not neces- 
sary for maintaining the stability of the dark matter; 
its purpose is to package the fine-tuning of the squared 
masses for the five additional doublets into a single 
tuning. This discrete symmetry also reduces the number 
of quartic couplings in the potential to that of the two 
Higgs doublet model. Any discrete symmetry group can 
be chosen, so long as it has a 5 representation (i.e., S§). 

The scalar potential for the six Higgs doublet model 

is 

V = -rr%\h\ 2 + TnZ\H i \ a + X 1 (\h\ 2 ) 2 + \ 2 (\H i \*) a 
+A 3 |/i| 2 |-ff,:| 2 + A 4 |/i t ff i | 2 + \ 5 ((fJHi) 2 + h.c.) 
+X 6 c ijk (hH}H j Hl + h.c), (1) 

where k — 1,...,5. Depending on the choice of 
discrete symmetry, there may be several couplings of 



the form |i?s| 4 ; one possibility is shown in the A2 term. 
The existence of the term proportional to A5 is necessary 
for the phcnomcnological viability of the model and 
forces the five-dimensional representation to be real. The 
term proportional to Xq is only allowed if the symmetry 
satisfies the following relation 

5<g>5 = l050---. (2) 

The couplings A5 and Xq lead to a physical phase and will 
induce CP violation in the self-interaction of the 5-plets 
after electroweak symmetry breaking. This does not 
alter the tree-level spectrum; because the self-coupling of 
the 5-plet is only affected at loop-level, the CP violation 
does not significantly alter the experimental signatures 
of the model. 

The field h acquires a vev and gives masses to the 
gauge bosons. In contrast, the field H5 does not acquire 
a vev and cannot have any Yukawa interactions with 
the Standard Model fermions. Expanding about the 
minimum of the potential, (h) — u/v2j with v = 246 
GeV. The Higgs 5-plets are 
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The physical masses of the particles at the minimum are 
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and must always be greater than zero. The lightest 
neutral particle, a°, serves as the dark matter candidate; 
in order that it not become the charged ± boson, the 
quartics must satsify 



A 4 - 2|A 5 | < 0. 



(5) 



Only the 5 option will be discussed here, but the results also 
hold for the 6 multiplet. 



The splitting between s° and a is proportional to A5 
and breaks the accidental U(l) symmetry. Results from 
direct detection experiments (see Sec. VAl require that 
the mass splitting between s° and a u be more than 
O(100 keV), which sets the limit 

|A 5 | > 10- 6 . (6) 

The experimental lower bound on the Higgs mass [23] 
constrains the value of Ai to be 

Ai>0.1. (7) 

Additional constraints on the quartics come from 
the requirement of vacuum stability. In order that 
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the potential (JlJ be bounded from below in all field 
directions, the couplings must satisfy 



Ai,A 2 > 

A 3 > -2^X 
A 3 + A 4 -2|A 5 | > -2^a7a 



(8) 



These conditions are for local stability of the potential 
at a given scale. If they are satisfied at all scales, then 
they correspond to absolute stability. 



III. RENORMALIZATION GROUP 
INFLUENCE ON LOW ENERGY SPECTRUM 

A. Gauge Unification 

Unification is the key motivation for introducing the 
five-plet H 5 and it is straightforward to check that the 
gauge couplings unify reasonably well with the addition 
of six or seven scalars to the SM (see [23] for two- 
loop RGEs) . In particular, when two-loop RGEs are 
evaluated, a threshold correction splitting a (fcrmionic 
and scalar) 5+5 by 7772/7773 ~ 30 is necessary to maintain 
unification for the six scalar case (where 777,2 and 1713 
are, respectively, the masses of a doublet and triplet). 
This is better than the case of seven scalars, where 
W2/7773 ~ 300. As a point of comparison, the threshold 
corrections for the MSSM require 7773/7772 — 20 [25] . 

The unification scale is given by 



*GUT = 2?T- 
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Mgut ^ 10 i4 GeV (9) 
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and the value of the gauge coupling at the GUT scale is 
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t GVT ~ 40. 



(10) 



If this theory is embedded in a simple SU(5) GUT, the 
resulting six-dimensional proton decay is 
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(11) 



which is far too fast. This implies that GUT- 
scale physics is non-minimal and must suppress gauge- 
mediated proton decay. Several approaches exist in 
the literature to deal with this task. One possibility, 
discussed in [2], is to embed the theory in a five- 
dimensional orbifold model. Proton decay is still al- 
lowed, but is highly suppressed due to the configuration 
of the fields in the extra dimensions. Trinification, a 
GUT based upon the group \SU(3)] 3 , provides another 
option because it completely forbids proton decay via 
gauge bosons [26] . 



B. Quartic Couplings 

The ability to discover the six Higgs dark matter 
candidate depends on its mass and its couplings to SM 
particles. The gauge interactions are fixed, but the 
couplings to the SM Higgs are model-dependent. These 
couplings must satisfy two requirements: perturbativity 
and vacuum stability. In addition, they must adhere 
to experimental constraints from Higgs and dark matter 
searches (see Sec. [H]). The model dependence comes into 
play when choosing the value of the quartics at the GUT 
scale. 

The most common understanding of how fine-tuning 
can give rise to Higgs and dark matter candidates near 
the electroweak scale invokes a landscape of vacua, each 
with its own values for couplings and masses. The string 
theory landscape allows for a great range of possibilities 
for the physical parameters of the theory and naturally 
leads to the question of what the typical values are 
in our neighborhood of vacua. The distribution of 
couplings is clearly a UV sensitive question and cannot 
be obtained by dimensional analysis because the quartics 
are dimcnsionlcss. Fortunately, there arc simple ansatze 
that lead to distinct weak-scale spectra. 

This section will explore two possible GUT-scale 
distributions of the quartics: parameter space democ- 
racy and susy. The couplings at the weak scale are 
obtained by application of the renormalization group 
equations. The resulting differences in weak-scale 
phenomenology for each of these distributions will be 
explored in Sec. [V] 



1. Parameter Space Democracy 

Perhaps the most obvious distribution of parameters 
is one where all couplings at the GUT scale arc equally 
probable - "parameter democracy." This measure favors 
large couplings of cither sign. If all the couplings are 
positive, they quickly run down to perturbative values 
and the initial boundary conditions of the quartics are 
not terribly important. When the quartic couplings start 
off negative, they can become asymptotically free and 
may have Landau poles. Furthermore, negative quartic 
couplings can lead to vacuum decay, especially when 
they have a large magnitude initially; thus, most of the 
parameter space in the negative direction is ruled out. 

Typically, the couplings approach a tracking so- 
lution rather rapidly [27] . Neither A5 nor Xq has a 
significant affect on the fixed point values of the other 
couplings. The gauge boson and top quark contributions 
also do not significantly affect the runnings in this region. 
With these observations, the beta functions may be 
approximated as 
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where 
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and Nh is the number of scalars added to the SM, 
in addition to the usual Higgs [28]. For the model 
considered here, Nh = 5. 

When Ai and A 2 are large at the GUT scale, the 
self-coupling terms in the beta functions ( 12 1 dominate 



and the low energy values for these couplings are approx- 
imately 
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Figure [I] (gray points) shows the weak-scale distribution 
of Ai and the effective coupling 
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(14) 



which parametrizes the interaction of the WIMP candi- 
date a® to the SM Higgs h°. The values of the quartics 
at the GUT scale were randomly sampled within the 
range: < Ai,A 2 ,|A 3 | < C(4tt), -1 < A 4 < 0, and 
I A5 1 < 2.1 They were then run down to the electroweak 
scale by applying the renormalization group equations. 
Despite the large range of possibilities at UV energies, 
the couplings are focused down to a narrow set at 
electroweak energies. Indeed, Ai and A 2 do not vary 



much from the values approximated in (13 1. The region 



of parameter space at the electroweak scale corresponds 
to 



0.1<Ai<0.3, 0<A 2 <0.1, 
-0.2 < A 3 < 0.4, -0.5 < A 4 < 0. 



(15) 



A5 renormalizes itself and thus remains small 
(IA5I < 0.1). With the assumption of parameter 
democracy, the model comes close to saturating the 
upper values for Ai and A 2 , having a small A 3 , and 
having a A4 that is close to saturating the lower bound. 
The acceptable range for the Higgs mass is 



114 GeV < 



< 200 GeV. 



(16) 



Higgs masses at the upper-end of this interval are 
preferred (see Fig. [TJ . 



t Quartics outside this range either give the same result for the 
low-energy spectra or, as in the case of A5, cause the couplings 
to run down to non-perturbative values. This range was chosen 
to maximize the sampling rate of the program. 
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FIG. 1: (Top) Distribution of Ai and A c ft = A3 + A4 — 
2 1 A5 1 at the electroweak scale obtained by solving the one- 
loop renormalization group equations for parameter space 
democracy (gray) and susy (black) boundary conditions at 
the GUT-scale. The couplings are focused down to a small 
range at the weak scale. (Bottom) The distribution of SM 
Higgs mass for the two sets of boundary conditions. The 
distribution of allowed Higgs masses is smaller in the case 
of susy boundary conditions as opposed to parameter space 
democracy conditions. 



2. Minimal Susy Boundary Conditions 

Another plausible set of boundary conditions are 
ones where supersymmetry is broken at the GUT scale 
and the dominant quartic couplings are those arising 
from D-terms. The simplest way of achieving the desired 
low-energy spectrum is if each low-energy Higgs doublet 
comes from a vector-like chiral superfield: $h and <f>£ 
for the Standard Model Higgs and 
five-plet of scalar dark matter. Specifically, 



and <&^ 5 for the 
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where [3 and (3$ are the orientation of the scalars inside 
the chiral superfields. The resulting D-term potential 
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has the following couplings 
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(17) 



A term must be added to the superpotential to generate 
A5 = 0. In order that this not alter the above 
relations significantly, it should be a small coupling. One 
possibility is to take the minimum value allowed by direct 
detection experiments, A5 ~ 1CT 6 . The gauge couplings 
at the GUT scale are 3q UT = 0.32, so the susy boundary 
conditions result in small couplings at the electroweak 
scale. In order to have neutral dark matter, A4 < 0, so 
cos2/3cos2/? 5 < 0. 

These couplings are a function of two angles and 
lead to a lighter Higgs and smaller mass splittings for 
the scalars than the case of parameter space democracy. 
This is apparent from Figure 1, where the black points 
show the allowed values of Ai and A g at the electroweak 
scale obtained using the susy boundary conditions. In 
this case, the Higgs mass falls within a much smaller 
range 



147 GeV < m h o < 159 GeV 



(18) 



and is lighter than the most probable Higgs mass for 
parameter space democracy. 



IV. DARK MATTER 



*5 ' 



A. Relic abundance 



The lightest neutral component of the H$ doublet, 
is a viable candidate for the observed dark matter 
and its mass may be estimated from standard relic 
abundance calculations. It is assumed that the a° is in 
thermal equilibrium during the early universe. When the 
annihilation rate of the a° is on the order of the Hubble 
constant, its number density 'freezes out,' resulting in 
the abundance seen today. 

When the freeze-out temperature is on the order of 
the mass splittings Am s o a o and Am^ n o, the presence of 
the additional scalars s° and 4»§ becomes relevant [29] . 
In this case, interactions involving the two other scalars 
as initial state particles are important in determining the 
relic abundance of a§ , which must fall within the WMAP 
region 0.099 < tldmh 2 < 0.113, where f^dm is the dark 
matter fraction of the critical density and h = 0.72 ±0.05 
is the Hubble constant in units of 100 km s _1 Mpc -1 ) 
30J . Typically, coannihilation has a significant effect on 
the allowed mass range of the relic. 



dn 

7ft 



The number density of a® is given by 
= -3Hn- ^ (<7ijVij)(niT 



»f ), (19) 



where cr^- is the sum of the annihilation cross sections of 
the new scalars Xi into Standard Model particles X 

Oij =^ i a{X i X i ^ XX). (20) 



The first term on the r.h.s. of equation (19) accounts 



for the decrease in the relic density due to the expansion 
of the universe; the second term results from dilution 
of the relic; from interactions with other particles. The 
annihilation rate depends on the number of scalars added 
to the theory in addition to the SM Higgs, Nh, through 
the interaction cross sections try. In general, Oi 



N h m a , so the mass of the dark matter scales as 



(21) 



Thus, in the non-resonance regime, the dark matter mass 
decreases with the number of electroweak doublets added 
to the SM. For this reason, the six Higgs doublet model 
gives lighter dark matter than the inert doublet model 

When m a o < 80 GeV, the only annihilation channel 
is to a pair of fermions. Because these cross sections 
tend to be rather small, H,^ m h 2 > 0.1. However, a 
resonance due to s-channel SM Higgs exchange causes 
a sharp decrease in the relic density ~ 80 GeV, bringing 
it within the WMAP experimental range. For m a o > 80 
GeV, diboson production is the dominant annihilation 
mechanism and keeps the abundance small. There is 
always another point in this large mass regime where 
^dm^ 2 ~ 0.1. Thus, the dark matter can take two 
possible mass values - one light (~ 80 GeV) and the 
other heavy (> 200 GeV). 

The relic abundance calculation was performed 
numerically by scanning over the parameter m 2 for each 
set of randomly selected quartic couplings. Figure [2] is 
a plot of the allowed mass of a® as a function of the 
SM Higgs mass. A broad range of values is allowed for 
the case of parameter space democracy, with m a o falling 
between ~ 200—700 GeV. For supersymmetric boundary 
conditions, the mass values range from ~ 200 — 400 GeV. 
An ~ 80 GeV dark matter particle is also allowed for 
both cases. 



B. Bounds from electroweak precision tests 

Electroweak precision tests place limits on the light 
mass range of the dark matter [31j . The contribution of 



5 



700 




g 200 



120 140 160 

SM Higgs mass (GeV) 



180 



200 



FIG. 2: Allowed mass of the LSP a% as a function of the 
Higgs mass for parameter space democracy (gray) and susy 
(black) boundary conditions. All points included in this plot 
fall within la of the electroweak precision data (Sect. IV B I 
and are consistent with LEP results (Sect. |V C| |. The SM 
Higgs can decay into a pair of WIMPs if m a o lies below the 
dashed red line. 



the new particles to the T parameter is given by 
N h 
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The expression for T(mi, 777,2) can be simplified if one 
assumes that the mass splitting Am = m 2 — mi satisfies 
Am/7711 <C 1- In this limit, 

/(Am) 4 ' 
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and the expression for AT reduces to 
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Because A5 is typically smaller than A4, AT is 
always positive. In the minimal SM, AT is driven 
more negative as the mass of the Higgs increases. 
The additional scalar doublet H§ compensates for this 
change, driving T positive. 

The S parameter also has contributions from the 
additional Higgs doublets [T!5] and is easily generalized 
to the case of six Higgses 
-1 
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When the mass splittings are small, 
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For the heavy dark matter candidate, the cor- 
rections to the S and T parameters fall well within 
the lcr electroweak precision data [T5]. Lighter dark 
matter can make significant contributions to the S and 
T parameters. Couplings that give rise to deviations 
in S and T that are more than lcr away from the 
measured values have not been used in the analysis of 
the experimental signatures of the model (Sect. |v| . 



V. EXPERIMENTAL SIGNATURES 
A. Direct detection 

Direct detection experiments provide a means for 
observing the dark matter relic when it scatters elasti- 
cally off atomic nuclei [52] . The WIMP can either couple 
to the spin of the nucleus or to its mass. The spin- 
independent contribution to the cross section usually 
dominates and bounds on its value are being set by ex- 
periments such as CDMS, DAMA, Edelweiss, ZEPLIN-I, 
and CRESST. 

In the six Higgs doublet model, there are two 
contributions to the spin-independent cross section. The 
first comes from an s-channel Higgs exchange described 
by the effective Lagrangian 
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Experimental results are usually reported in terms of the 
cross section per nucleon, which in this case is 
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Because 
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This cross section scales as Nh (see Eq. 
o n oc m~ 2 , the lighter dark matter candidate will have a 
stronger signal than its heavier counterpart. Figure [3j 
shows the cross section per nucleon for the case of 
parameter space democracy (gray) and susy (black) 
boundary conditions. The current CDMS II run is 
sensitive to the lightest WIMPs predicted by the model. 
A larger portion of the parameter space is within the 
testable reach of the proposed SuperCDMS experiment 
[55] . The lower dashed line on the plot is the expected 
limit from Phase C of SuperCDMS. 

Another contribution to the spin-independent cross 
section comes from the inelastic vector-like interaction 



G 



10" 



10" 



10" 





















( 


3DMS II 2007 


















1 








SuperCDMS 














• 






: 


















! 











200 



400 600 
WIMP mass (GeV) 



800 



1000 



FIG. 3: Cross section per nucleon for the case of parameter 
space democracy (gray) and susy (black) boundary condi- 
tions. The lightest WIMP candidates will be tested for at the 
current CDMS II run (upper dashed line). The third phase of 
SuperCDMS (lower dashed line) will probe a greater region 
of the parameter space. 



a + p — > s° + p, which is mediated by an off-shell Z°- 
boson. In general, such inelastic transitions provide a 
means to reconcile DAMA's detection of relic-nucleon 
scattering, which conflicts with CDMS's null result [M] , 
Consistency with the experimental results requires a 
mass splitting Am s o o ~ 100 keV between the two 
lightest scalars [33]. 



B. Indirect detection 

A concentration of WIMPs in the galactic halo 
increases the probability that they will annihilate to 
produce high-energy gamma rays and positrons |36j . 
The gamma ray signal is of particular interest because 
it is not scattered by the intergalactic medium; thus, 
it should be possible to extract information about the 
WIMP mass from the spectrum. 

Monochromatic photons can be produced when the 
WIMP annihilates to produce 77 and Z°-f. The dom- 
inant mechanisms that contribute to this annihilation 
depend on the DM mass regime. The light dark matter, 
for example, annihilates primarily through s-channcl 
Higgs exchange with a one-loop h°^X vertex (X = 7, 
Z°). The main contributions to the loop come from the 
W ± boson, the top quark, and the <j) 5 five-plet. Other 
box diagrams are suppressed. The WIMPs are highly 
non-relativistic and their annihilation cross section in the 
light mass regime is nearly 
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where u is the relative velocity between the initial two 
WIMPs and s sa 4m^ - The general expressions for the 



decay widths of the Higgs boson into a 77 and jZ° final 
state are found in (37J |38] . 

The case of the heavy dark matter is significantly 
different |39j . In this regime, the dominant contribution 
comes from the box diagram with three (frg and one 
W + in the loop. When the a 5 and cf>^ are nearly 
degenerate and m a o 3> my/± , there is an effective long- 
range Yukawa force between the </>5~</>5 pair in the loop 
that is mediated by the gauge boson: 



V(r) 



-0L2- 



(31) 



As a result, the pair of charged scalars form a 
bound-state solution to the non-relativistic Schrodinger 
equation. The optical theorem is used to obtain the s- 
wave production cross section for the bound state: 



a(a 5 a 5 -> 0jj> 5 )u ■ 



N h m 2 




2mAr 



(32) 

Multiplying this by the decay width of the bound state 
to two photons (or, jZ°), gives the total annihilation 
cross section 



a(a 5 a 5 -> 77)?/ 



2ira 2 a% 
N h m 2 



w± 




2toAto 



<£±a° 



(33) 

This cross section does not depend on m a o (to zero- 
th order in the mass splittings) and, as a result, is 
significantly enhanced in the heavy DM mass region. 
This enhancement is critical; because of it, the heavy 
mass DM may be visible in gamma ray experiments. 
Additionally, the only parameter dependence comes in 
through the mass-splittings, which are small. Therefore, 
there is not much spread in the range of allowed cross 
sections. 

The monochromatic flux due to the gamma ray final 
states observed by a telescope with a field of view Afi 
and line of sight parametrized by W = (0, <j)) is given by 



$ = a 



( VyXi 



where 



V 1 P b 



y 77 



)( 



100 GeV\ 2 



J(*,Afi)Afi, (34) 



i-9„ m -2„-l 



C lZ a = 5.5 x 10~ xo cm~ 2 s _1 



(35) 



and the function J includes the information about the 
dark matter distribution in the halo. Note that the flux 
is independent of Nh- For the NFW profile, J ~ 10 3 for 
Af2 = 10 -3 [35]. Other profile models exist with either 
more mildly /strongly cusped profiles at the galactic 
center [37 . Depending on which model is chosen, J can 
be as small as 10 or as large as 10 5 for Af2 = 10~ 3 . 
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In this work, the moderate NFW profile will be used, 
however the result can easily be scaled by two orders of 
magnitude to get the predictions for other halo profiles. 

The expected monochromatic flux for the 77 line 
is shown in Figure [4] (assuming JA£l = 1). The 
estimated flux is right beneath the sensitivity limits of 
the ground-based HESS detector (green line) and space- 
based GLAST telescope (red line) [371110]. The results 
for the 7Z line are similar for low DM masses and are 
enhanced by an order of magnitude for masses greater 
than 200 GeV, putting it within the reach of HESS. 
Given the two order of magnitude uncertainty in the flux 
coming from the details of the halo profile, the gamma 
ray line is an interesting signal for both current and 
upcoming experiments. 



C. Collider Signatures 

It is possible that the scalars of the Higgs 5-plet 
were produced at the e + e~ collider LEP with ^/s ~ 200 
GeV via the processes 



^5 



and 



(36) 



LEP placed limits on the production cross sections for 
the neutralino and chargino |41j and these bounds can 
be directly translated to the processes in (36 1. By doing 



so, approximate limits on the masses of the scalars 05, S5, 
and 4> 5 can be deduced. The charged Higgs 5-plet cj>f 
is ruled out for masses below ~ 90 GeV and the neutral 
scalar is ruled out for masses between ~ 100 — 120 
GeV (depending on the mass of 05). Fig. [5] summarizes 
the important constraints on the mass splittings A 
and Am s o a o for the light dark matter. 



The heavy dark matter candidate (m a o > 200 
GeV) could not be produced at LEP. In addition, 
its contributions to the electroweak parameters always 
fall within the ler experimental bounds. Thus, the 
main constraints on the mass splittings in this region 
of parameter space come from vacuum stability and 
perturbativity. Typically, A Am s o a o < 20 GeV 

in the heavy dark matter regime. 

One of the most promising discovery channels for 
the 5-plet scalars at the Tevatron and LHC is the width 
of the SM Higgs. The Higgs can decay into 05,55, 
or <fif , in addition to the SM modes. In Fig. V2\ all 
points below the red dotted line satisfy m^a > 2m a o; 
here, the SM Higgs can decay into the dark matter. 
This decay channel is open for significant portions of 
both the parameter space democracy and susy boundary 
condition cases. For small mass splittings, decays to S5 
and 4> 5 are also possible, though they are subdominant. 

The contribution of the new invisible decays to the 
width of the Higgs is 



r. — 

x inv 



N h v 2 
327rTO/,o 




Fig.[6]plots ri nv as a function of the SM Higgs mass. 
The (top) line is the width due to the SM decay modes 
|4"2"] . For points above the line, the invisible decays 
into the 5-plet scalars are the dominant contribution. 
However, even when rj nv ~ O.IFsm, it should be possible 
to detect the additional decay modes at the Tevatron or 
the LHC. 




100 200 300 400 500 600 700 
WIMP mass (GeV) 




20 40 60 

Am sa =m s -jn a (GeV) 



FIG. 4: Approximate flux from dark matter annihilation 
in the galactic halo via a® a® — > 77 for parameter space 
democracy (gray) and susy (black) boundary conditions. The 
dashed lines indicate the sensitivity of GLAST (red) and the 
ground-based detector HESS (green). The NFW profile was 
used. 



FIG. 5: Mass splittings for the light dark matter (m a o ~ 80 
GeV). LEP excludes the region of intermediate S5 mass. The 
region on the upper left is excluded by electroweak precision 
results, while that on the lower right is excluded by vacuum 
stability. Results are shown for parameter space democracy 
(gray) and susy (black) boundary conditions. 



8 



The scalars may also be produced at the Large 
Hadron Collider via interactions like 



VV - 


- agag - 


> Z* + 




VV - 








vv - 


- sUf ~ 


-4 Z*W* +fi T 




vv - 




-> W* + fi T . 


(38) 



The vector bosons are always off-shell because the scalar 
mass splittings are less than 80 GeV (see Fig. |5| and, 
after using their leptonic branching fraction, it will be 
challenging to detect this signal in the presence of a large 
background. 

As an example, consider the first two processes in 
( 38 1 , which both result in opposite-sign leptons plus $t 
after the decay of the gauge bosons. The production 
cross section <7 pro d for s^a^ and at the LHC was 

calculated using MadGraph [35] for a sample point in 
parameter space and is plotted as a function of WIMP 
mass in Fig. [7] The cross section for the SM background 
(thick line) is 



^background 



a(pp- 

+<j(vp 



WW)Bt(W 
-> ZZ)Bt(Z - 



i+r 



)Br(Z 



vv). 
(39) 



The signal cross section may be estimated as 

^signal ~ Br(Z — > l + r)cr plod , (40) 

where the branching fraction is about 1% (dashed line). 
The ratio of signal to background is about 1:10 for the 
low mass dark matter. At higher mass, it is about 1:1000. 
This estimate indicates that it may be possible to see the 
signal for the low-mass DM region, if appropriate cuts 
are placed. 



VI. CONCLUSIONS 

In this paper, a minimal extension of the Standard 
Model was presented that lead to gauge unification and 
a dark matter candidate. An electroweak doublet Hi in 
a 5 of a global discrete symmetry was introduced. One 
of the new five-plet of particles is light and neutral and 
serves as a good dark matter candidate. The addition of 
six scalars to the SM leads to gauge coupling unification 
and fixes the number of electroweak doublets. 

The six Higgs doublet model has distinct signatures 
for direct detection, indirect detection, and collider 
experiments. Typically, the light mass range (m a o ~ 
80 GeV) has the most promising signals, and will be 
tested for by GLAST and CDMS. In addition, it can be 
produced by decays of the SM Higgs at the Tevatron or 
LHC. The heavier candidates (m a a > 200 GeV) are more 
difficult to see, but lie within the sensitivity of the HESS 
gamma ray detector and the next-generation direct 
detection experiment, SuperCDMS. Direct production 
of these heavier candidates at colliders is challenging 
due to large Standard Model backgrounds from di-boson 
production, though further study is needed to determine 
whether appropriate cuts can reduce these backgrounds. 

Throughout this discussion, it has been assumed 
that an exact discrete symmetry exists to keep the full 
five-plet of dark matter light under one fine-tuning. 
Discrete symmetries can arise in string theoretic con- 
structions (e.g., see [H]), and the existence of these sym- 
metries is critical for viability of this particular model. 
If this requisite symmetry is relatively common, then a 
single fine-tuning of the scalar mass is comparable to 
the tuning necessary in the "minimal model" described 
in [21 [TS] . In general, this class of minimal models is 
not as economical in terms of fine-tuning as split susy, 
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rod = "-(pp->a s s 5 ) 


+ o-(pp->f 5 V 5 ) 


















- ^ ^"signal 















200 300 400 

WIMP mass (GeV) 



FIG. 6: Width of the SM Higgs decay into the H5 scalars 
for the case of parameter space democracy (gray) and susy 
(black) boundary conditions. The top line shows the width 
of the SM decay modes, Fsm- The bottom line is O.IFsm- 



FIG. 7: LHC production cross section <T pro d for agss and 
05 05 , assuming Am s o a o = 10 GeV, Am^± a o = 15 GeV, 
m h o — 120 GeV, and A3 = 0.3. The dotted line is the signal 
cross section. The thick black line is the cross section for the 
SM background (see text). 
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where the desired mass spectrum is obtained by having 
the R-symmetry breaking scale be small. However, these 
minimal models give rise to interesting phenomenology 
that will be tested in upcoming experiments. 
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